On the reliable measurement of the gluon condensate in lattice QCD 
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We propose to calculate the gluon condensate in lattice QCD in an indirect way by extracting it from the 
correlator of hadronic currents of heavy quarks. Moments (derivatives with respect to momentum at vanishing 
momentum) of the vector and pseudoscalar correlators are evaluated. The contribution of the continuum spectrum 
in addition to the low-lying resonance is observed in accordance with qualitative expectations. Finite size effects 
are observed and shown to be insignificant for heavy quarks and moderate lattices. A practical definition of the 
nonperturbatively renormalized heavy quark mass is given. 
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Introduction 

Extensive analysis of the two-point correlators 
of hadronic currents by means of the QCD sum 
rules indicates that vacuum expectation values 
of the quark and gluon fields receive large non- 
perturbative contributions (condensates) [1]. The 
existence of condensates implies a particular form 
of the crossover between the domain of asymp- 
totic freedom and the domain of large distances, 
dominated by the low-lying resonance. There is 
certain range of intermediate virtualities where 
the perturbative series in powers of a s is still un- 
der control, whereas nonperturbative effects show 
up, represented, in some channels, by the first 
few terms of the operator product expansion. It 
was recently suggested [2] that lattice simulations 
must be useful to study those intermediate dis- 
tances in the physics of light quarks. We find it 
very advantageous to study deviations from the 
asymptotic freedom in the case of correlators of 
hadronic currents of heavy quarks. It is in this 
case that the gluon condensate was first intro- 
duced [1]. Such a way is an alternative to the 
direct evaluation of the gluon condensate as an 
average of the local operator, initiated in [3]. 

Consider the polarization operator of the elec- 
tromagnetic currents of charmed quarks: 



IV = i / dxe'«* < 0\T(j^x)jA0)) |0 > 



where j^ = cj^c. The function Ii-(q 2 ) is com- 
putable in perturbative QCD at the origin q 2 = 
since that point is far away from the leading 
threshold due to a pair of charmed quark and 
antiquark. Ii-(q 2 = 0) is a function of the renor- 
malized heavy quark mass m c (/u) and a s (ji). The 
relevant renormalization point here /j, = 2m c is 
high enough to ensure the validity of the QCD 
perturbation theory: a s (2m c ) ~ 0.2. Define mo- 
ments of the polarization operator (1) as: 
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Starting with n = 2 the moments Mn§ are phys- 
ical quantities: all the ultraviolet singularities 
are adsorbed by the bare Lagrangian parameters. 
Higher moments n >> 1 correspond to smaller 
virtualities, which leads to deviations from the 
perturbative series before that perturbative se- 
ries grows out of control [1]. For n = 4 -=- 7 
those deviations are well described by the few 
terms of the operator product expansion, start- 
ing with the quadratic gluon condensate G 1 - 2 - 1 = 
< 0|(a s /7r) G^ V G^ V |0 >. For the ratios of neigh- 
boring moments r n = M n j ' M n -\ one has on 
the two-loop level : 
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where {a n , b n ,c n } are known numbers. The point 
we are making is that the l.h.s. of (3) is di- 



rectly accessible to the Monte-Carlo lattice sim- 
ulations, and since the perturbative series in the 
r.h.s. of (3) is reliable the nonperturbative term 
represented by the gluon condensate is to be ex- 
tracted in a reliable way. Lower moments (7*34) 
correspond to short distances and are dominated 
by the first two terms of (3), where a n corre- 
sponds to the one-loop contribution (asymptotic 
freedom term) and b n represents a small two- 
loop correction. Since the one-loop diagram is 
a function of only one parameter - the renormal- 
ized heavy quark mass, that parameter can be ex- 
tracted from the first moments. In this way, using 
Monte-Carlo data on a finite lattice one can deter- 
mine the nonperturbatively renormalized heavy 
quark mass in terms of the gauge invariant cor- 
relator without any reference to the pole mass. 
This mass is to be used later to extract the last 
term in the r.h.s. of (3) from the higher moments. 
For any given Monte-Carlo evaluation of the l.h.s. 
of (3) the perturbative part of the r.h.s. of that 
equation is to be built in terms of renormalised 
parameters of the improved lattice perturbation 
theory [4]. 

The contribution of the gluon condensate in (3) 
is about 10% relative to the leading first term 
~ a n in the charmonium case. However, we are 
not forced to consider the phenomenological value 
of m c . Adjusting the hopping parameter k so 
as to decrease m c from l.SGeV to, say 0.9GeV 
where perturbative series in a s (0.9GeV) is still 
under control [1], one can easily increase the third 
term in (3) up to 20%. The efficiency of this trick 
is due to the fact that the perturbative terms of 
the r.h.s. of (3) depend on m c logarithmically, 
while the gluon condensate contribution is rapidly 
varying power correction mj 4 . 

On the lattice one has to deal with the orig- 
inal correlator Jl^ v {q 2 ) rather than T[(q 2 ). We 
therefore introduce moments M n as derivatives 
of n /U/U (g 2 ) rather than II((/ 2 ) in (2). The new 
moments are related to the initial ones (2) in a 
simple way M n = — 3.M n -i, so the ratio r n for 
the new moments will be f n+ \ = r n . To know 
the Fourier transform of the correlator (1) from 
the lattice simulations to the accuracy needed for 
the detection of the power corrections does not 



seem to be feasible. For this reason we evaluate 
not the Fourier transform Jl^ v {q 2 ) at nonzero q 2 , 
but rather the moments themself directly in the 
;c-space. At q 2 = one can get rather simple ex- 
pressions for the moments in terms of the original 
correlator IL^^x): 
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Monte- Carlo data 



We used the clover- and tadpole-improved 
quenched Wilson fermionic propagators of [5], 
obtained on the 8 3 x 16 lattice blocked from 
32 3 x 64 S77(3) configurations at j3 = 6. We 
have calculated the correlator (1) and the ra- 
tios of the moments (3) for vector (Fig.l) and 
pseudoscalar (Fig. 2) currents (j = icy$c in (1)) 
for three values of the hopping parameter k = 
(0.9830,0.1031,0.1111). Statistical errors are no- 
ticeable but small for already 6 propagators (n = 
0.1111 - the upper curves) and negligible for 10 
propagators (n = 0.0983 - the lower curves). The 
given values of k correspond to not heavy quarks, 
so our data for r n should be compared to the QCD 
perturbation theory as in (3). We rather analyse 
the form of the spectral density p(s) of the corre- 
lator (1) by looking at the functions r(n). 

The simple and phenomenologically very suc- 
cessful ansatz for p(s): 

p(s) = fm 2 e J(s - m 2 es ) + —0{s-s o ) (5) 

47T^ 

with a single low-lying resonance of mass m res , 
residue / and perturbative continuum with some 
effective threshold s > m 2 es gives the ratios 
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At large n the ratios r n reach a constant l/m 2 es , 
which is the resonance contribution. That con- 
stant is approached from below due to the con- 
tinuum contribution, which is just what one can 
see on Fig.l. As n grows the ratios try to reach 
the asymptotic constant value from below, but 



at some point fail. That "crash" -point is a fi- 
nite size effect. One cannot put arbitrarily high 
moments on a finite lattice because higher mo- 
ments correspond to lower virtualities, larger dis- 
tances. At some point they become too sensitive 
to the boundaries of the box. This interpreta- 
tion is supported by the fact the "crash" -point is 
clearly K-dependent. At greater k, smaller m res 
(upper curves) the "crash" happens earlier. Figs. 
1,2 demonstrate that one can easily accommodate 
heavy fermions (n < .0983) on typically simulated 
lattices as far as the ratios of interest r^ -r- J*7 [1] 
are concerned. 
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continuum threshold s is significantly greater 
than resonance mass in that channel. For heavy 
quarks the continuum contribution will be more 
visible for small n and allow for the comparison 
with the predictions of perturbative QCD. The 
Monte-Caro data with smaller lattice-spacings 
and smaller k are necessary to see the slope of 
the continuum contribution with good resolution 
at small n. 
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Figure 2. Same as in Fig. 1 in the pseudoscalar 
channel. 
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Figure 1. Ratios r n for three different hopping 
parameters in the vector channel. 



The observed "asymptotic" values of the ratios 
do follow the law ~ l/m^ es . For instance, in the 
vector channel the height of the upper curve is 
1.69 greater than the height of the lower curve. 
Whereas for the corresponding masses one has [5]: 



1.68 



(7) 



m 2 res (K = 0.0983) 
ml es (n = 0.1111) 

In the pseudoscalar channel (Fig. 2), where the 
resonance (pion) mass is smaller, the ratios r n are 
flat near the origin, which means that the contin- 
uum is not seen. That is because the effective 
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